Nearby quasar remnants and ultra-high energy cosmic rays 
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Abstract 

As recently suggested, nearby quasar remnants are plausible sites of black-hole based compact dy- 
namos that could be capable of accelerating ultra- high energy cosmic rays (UHECRs) . In such a model, 
UHECRs would originate at the nuclei of nearby dead quasars, those in which the putative underlying 
supermassive black holes are suitably spun-up. Based on galactic optical luminosity, morphological type, 
and redshift, we have compiled a small sample of nearby objects selected to be highly luminous, bulge- 
dominated galaxies, likely quasar remnants. The sky coordinates of these galaxies were then correlated 
with the arrival directions of cosmic rays detected at energies > 40 EeV. An apparently significant cor- 
relation appears in our data. This correlation appears at closer angular scales than those expected when 
taking into account the deflection caused by typically assumed IGM or galactic magnetic fields over a 
charged particle trajectory. Possible scenarios producing this effect are discussed, as is the astrophysics 
of the quasar remnant candidates. We suggest that quasar remnants be also taken into account in the 
forthcoming detailed search for correlations using data from the Auger Observatory. 



PACS number(s): 98.70.Sa, 98.54.-h 



^ 1 Introduction 

• i— i . 

^ Different experiments over the past few decades have detected several giant air showers, confirming the 
arrival of cosmic rays (CRs) with energies up to a few hundred EeV (1 EeV = 10 18 eV)[jl]]. The nature and 



origin of these energetic particles remain a mystery |2[. CRs, however, can not travel unaffected through 
inter-galactic space. The thermal photon background becomes highly blueshifted for ultra-relativistic pro- 
tons, and the reaction ^7 — > A + — > tt°p, and similar others, effectively degrade the primary proton energy. 
This provides a strong constraint on the proximity of CR-sources, discovered early on by Greisen, and by 
Zatsepin & Kuz'min, and referred to as the GZK cutoff |J. Specifically, fewer than 20% of 300 EeV (100 
EeV) protons can survive a trip of 18 (60) Mpc. For nuclei and photons the situation is, in general, more 
drastic f|. 

The possible solutions to this puzzle appear to fall into three broad categories, viz: 1) There are many 
nearby sources (e.g., based on a substantial present-epoch population of galaxies hosting core supermassive 
black holes). 2) There are only a few nearby sources (e.g. ||]) and particles are isotropized by strong de- 
flections in Galactic and/or extragalactic magnetic fields of micro-Gauss strength, close to existing upper 
limits (e.g., see ||], also |7|). 3) The particles somehow travel relatively unhindered through the Cosmic 
Microwave Background radiation field, either by virtue of being some exotic new kind of weakly photon- 
interacting entity, or by violating the Lorentz symmetry of special relativity (e.g. S), or by an as yet 
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unknown other effect. 

The arrival directions of the primary particles could be a useful source of information about the ori- 
gin(s) of CRs. However, due to the small number of events, particularly at the highest energies, correlation 
studies are to be considered preliminary. For instance, the first five events observed with E > 80 EeV did 
in fact point toward high redshift radio-loud quasars, astrophysical environments that could well accelerate 
CRs above the GZK energies via shock mechanisms (see Farrar and Biermann, in Ref. ||). However, with 
the inclusion of subsequent data, this association now seems to have disappeared []|]. 

In this communication we explore whether UHECRs can be in any way related to present-epoch su- 
permassive black holes. To that end, we analyze what we can learn from the arrival directions of cosmic 
rays concerning their possible correlation with nearby massive dark objects, candidate QRs (quasar rem- 
nants), and their underlying astrophysics .Q The compact dynamo model has been proposed as a natural 
mechanism for accelerating cosmic rays in such environments fllOfl . In this model UHECRs are produced in 
nearby dead quasars harboring spinning supermassive black holes. The required emf is generated by the 
black hole induced rotation of externally supplied magnetic field lines threading the horizon. The observed 
flux of CRs would apparently drain only a negligible amount of energy from the black hole dynamo, and 
particles up to at least 100 EeV are expected. It is then interesting to ask if we are able to see any 
correlation between the sky position of these QRs and those of the highest energy CRs, assuming different 
simple configurations for the intervening magnetic field. 



2 Quasar remnant candidates 

The Nearby Optical Galaxy (NOG) catalog of Giuricin et al. [O] is a complete magnitude-limited (cor- 
rected blue total magnitude B < 14), distance-limited (redshift z < 0.02) sample of several thousand 



galaxies of latitude \b\ > 20°, with their morphology T also provided [14|. We shall impose very re- 



strictive selection criteria, those necessary for obtaining candidate objects providing the most favorable 



setting for a black hole based compact dynamo model of UHECR production [ 10 1 . The key goal of our 
strategy is to select a priori an NOG subsample of galaxies that are likely to be quasar remnants. There- 
fore, we are seeking optically bright galaxies whose luminosities are bulge dominated (e.g., giant ellipti- 
cals). Imposing a GZK-related horizon of ~ 50 Mpc assures that there are no quasars in our sample. 
We already know that the best determined (and most massive) black hole nuclei tend to be associated 
with bulge luminosities corresponding to Mb ~ —21 and brighter (for tables, see Kormendy's website at 



http: / /chandra.as.utexas.edu/^kormendy/wwwbhtable-tech or the second reference in [27]). 



Firstly, then, we shall impose a cutoff in redshift, requiring all galaxies in our sample to be within 
z < 0.01, which for Hq = 75 km/s/Mpc corresponds to 40 Mpc. Since the CRs detected are mostly in the 
northern hemisphere (except those coming from SUGAR, which are not used in the present analysis), we 
shall impose the restriction that all QRs have equatorial latitudes north of —10 deg. We further require 
the absolute blue magnitude to be brighter than Mb = —21 (for Hq = 75 km/s/Mpc) and an RC3 (Third 
Reference Catalog of Bright Galaxies) morphological type T < —3. A more negative T indicates greater 
bulge prominence. In order to correlate the UHECR arrival directions with their putative QR origins, we 
require that these candidate QRs not lie within rich clusters (those having more than 50 members). Some 
of the most massive QRs could reside in rich clusters of galaxies (e.g., all four QRs considered by Boldt and 
Loewenstein in Ref. [[n]] are in rich clusters), but the magnetic field strength in those clusters is presum- 

lr The term "quasar remnants" was introduced by Chokshi and Turner JH]] to describe the present-epoch population of dead 
quasars harboring supermassive black hole nuclei. 
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Table 1: Sample of QR candidates, columns are the B1950 coordinates (a, S), name of the galaxy, richness 
of the group to which the galaxy pertains (number of members), galactic longitude and latitude, mean 
redshift, corrected redshift (against Hubble distortions, using models of the peculiar velocity fields Marinoni 
et al. |l5[), morphological type code (RC3), and absolute magnitude (for Hq =75). This table is based on 
data from the NOG catalog of Giuricin et al. provided to us by C. Marinoni [|i~2f . 
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ably several micro-Gauss [13]; hence, UHECRs originating in such QRs would be extremely deflected away 
from their sources. In addition, there are only a handful of galaxies within large clusters in the Giuricin et 
al.'s sample that survive all other constraints in order to be declared plausible quasar remnant candidates. 
Using these selection criteria, we obtain a sample of 12 candidate QRs. These candidate QRs are listed in 
Table 1, and their astrophysical properties are analyzed below. 

The UHECR sample used is that obtained with AGASA above 40 EeV |jlql. There are 38 such events 
at \b\ > 20°; the angular precision of their arrival directions was estimated as a circle of radius 1.6 degrees. 
For energies > 100 EeV we also consider the 7 events observed at \b\ > 20° compiled by Sigl et al. ||, 
4 of which were obtained with AGASAj^J We note that the statistical test we shall report was blind, i.e. 
we did not know beforehand if any of the QR galaxies were coincident with high energy cosmic rays. The 
simulation technique is the same as that used in Refs. g, [1], |l[. 



Looking for superposed QRs in the cosmic ray error circles, we find the results reported in Table [2|. 
When we consider the AGASA sample with E > 40 EeV, there is an excess in the real result which ap- 
proaches 3a. This could argue in favor of a correlation between our sample of QRs and UHECRs at small 
angular scales. This is not the case for the CRs with energies above 100 EeV, for which we found only 
an apparent excess for the nominal error circles, not significant enough to give to it any confidence. A 
summary of these results is given in Figure [l]. There we show the number of real coincidences in rings 
surrounding the position of the central QR as a function of the internal radius of the ring, A. Each ring 
has an angular width of 1.6 deg, and the number of coincidences is shown with angular offsets between A 
and A + 1.6. An excess of coincidences is clearly evident for the innermost radii. At larger radii, the real 
superpositions are compatible with the random coincidences. This result stands disregarding the value 
chosen for the size of the ring. 



2 See Table 1 in Q for this compilation and the associated error estimates. 
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Figure 1: Real coincidences in rings surrounding the position of the central QR as a function of the internal 
radius. In black we show the random result expectation, whereas in light colour we show the real result, 
x-axis is given in degrees. Error bars are la. 

Note that when increasing the error circles, no additional QRs are involved in the coincidences, but 
instead, there are new CRs involved. This is what makes the real result increase. Details of the coinciding 
pairs are given in Table ||[ It is risky, of course, to evaluate the level of significance that we should attribute 
to this correlation, especially when so few events/ objects constitute the samples involved. The 12 galaxies, 
however, are selected a priori on a restrictive physical basis, aimed at identifying QRs harboring super- 
massive black holes. Furthermore, we can see that the correlated pairs identify an apparently preferred 
region of the sky, as Figure § shows. The expected number of random coincidences decreases from 0.40 to 
0.14 when the region of the sky is restricted to an area given roughly by 160 < a < 200 and 50 < 5 < 80. 
This is ~ 5a below the real result of 2 coincidences. This latter result, however, is obtained a posteriori, 
and its importance is thereby reduced. 

If we consider that charged particles (protons) are accelerated in these QRs and then travel through the 
intergalactic and Galactic magnetic fields towards the Earth, we expect their trajectories to be deflected. 
When the Larmor radius of a particle (rx — 10 2 Mpc E20/B-9) is much larger than the coherence length 
of the magnetic field £ C oh, the characteristic deflection angle 8 from the direction of the source, located at 
a distance D, can be estimated assuming that the particle makes a random walk in the magnetic field (see 

G^f fef (£)• 

where E20 is the energy of the particle in units of 10 20 eV, and U_g is the magnetic field in units of 10 -9 G. 
We can see that scattering in large scale magnetic irreg ularities O (nG) @ are enough to bend the orbits 
of trans-GZK protons by about 4 deg in a 50 Mpc traversal. As exhibited in Table 3, the CR angular 
offsets observed for these QRs are much smaller than 8. For a variety of assumed magnetic field scenarios, 
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Table 2: Coincidences with the 34 AG AS A CRs with energies of E > 40 EeV and the given constraints in 
both 5 and b. The 12 QRs given in Table 1 are considered. The second panel shows similar results, but in 
this case, we consider the coincidences with 7 UHECR with E > 100 EeV in the same latitude range. A 
real result equal to 2 means that there are 2 different cosmic rays coinciding with the quoted QRs. 
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Table 3: Superposed pairs of QRs and CRs. Nominal errors are considered. 9(E) is given for the nominal 
values in its definition. The offset angles for NGC 3610 and NGC 3613 refer to the same CR, which lies 
within both error circles. 
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6 is often substantially larger than the estimated AGASA measurement error of 1.6 degrees. If, in fact, 
due to the value of the inter-galactic magnetic field, such is the case, it would then indicate that 

1. the number of apparent QR/UHECR associations is no longer above random expectations (see Figure 
1 for large value of A) , and 

2. the a priori probability is relatively very small for offset angles as little as those actually observed 
(Table 3). 

Should the apparent clustering of correlated pairs be supported by future data, what are the viable sce- 
narios under which this could occur? One possibility is to consider that the intergalactic medium between 
Earth and the three apparently 'clustered' QRs is sufficiently different from the intergalactic medium in 
front of the remaining nine objects that are much more uniformly distributed on the accessible sky. For the 
deflection of an energetic (60 EeV, the mean of the two CR energies in Table 3) proton in traversing 34 Mpc 
(the mean of the 3 QR distances in Table 3) to be less than a degree, we would need B < 2 x 10~ 10 
G, which appears to be not as drastic a difference from the canonical nano-Gauss B field and coherence 
length of 1 Mpc that are usually assumed. Indeed, the IGM B field is likely to be an order of magnitude 
less than a nano-Gauss in voids comparable in size to the GZK horizon (Peter Biermann; 2002, personal 
communication). Independent estimates of the Intergalactic Medium (IGM) magnetic field towards these 
galaxies would prove very useful in assessing this explanation. 

It is important to note that, in some directions, the magnetic field of our own galaxy could well lead to 
a deflection of up to several degrees for primaries with energies below 60 EeV; see, for instance, Table 1 of 
Ref. p2|| . A recent study of this issue was presented by Alvarez-Muhiz et al. ]23|. However, it is important 
to realize that the possible filamentary topology of the Galaxy's magnetic field (Gerrit Verschuur; 2002, 
personal communication) would likely allow some directional windows, albeit narrow, where the deflection 



5 




Figure 2: Sky distribution (galactic coordinates) of the selected QRs. The lighter dots stand for those 
QRs which exhibit superposition with UHECRs. Two closely neighboring QRs (in projection) are actually 
almost indistinguishable in this plots (NGC 3610 and NGC 3613). 



of an UHECR could be much less than typical. Deflections due to the magnetic fields would of course be 
avoided if the primary were a photon, generated in the neighborhood of the QR via an accelerated charged 
particle interaction. A photon primary is ruled out for the highest cosmic ray energy event detected with 
Fly's Eye ||], but there is yet some discussion of the likelihood of having some photons as primaries in the 
AGASA sample [24|. Protheroe and Johnson (see second reference in Q) have studied the propagation of 
ultra-high energy gamma rays in the intergalactic magnetic field in detail. They have shown that if the 
intergalactic magnetic field is low enough, say B ~ 1(T 10 G, a photon can survive a distance comparable 
to those of the QRs in our sample and be a plausible primary for the observed showers. A new neutral 
hadron could obviously be another extreme possibility. However, as we have seen, the involved QRs are all 
within a relatively small region of the sky; their angular separation is clearly less than the average angular 
separation between neighbors in the remaining sample of nine. Why would sources of neutral hadronic 
CRs be restricted to a limited angular region of the sky? Albeit suggestive, we warn the reader that the 
effect observed here might be no more than an artifact resulting from the small number statistics available 
at this time. It is, however, interesting to study the underlying astrophysics of our candidates as possible 
cosmic rays emitters. 



3 Underlying Astrophysics 

Table 4 exhibits the physical attributes of our sample of QR candidates, viz: morphological type, mass 
of supermassive black hole nucleus (Mbh)j galactic bulge mass (Mgb)> black hole accretion rate m (in 
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Eddington units), the maximum resulting proton energy (E maK ) expected, and lower limits to the proton's 
radiation length (A) for losses arising from photo-pion production in ultra-relativistic collisions with the 
ambient electromagnetic radiation (photon) field. That the nuclei of these galaxies are indeed dark as 
well as massive may be appreciated by noting that the radio luminosity {yL v at 6 cm) within a 5 arcsec 
beam ]25| is in no instance greater than 10 -8 (in Eddington units); the ROSAT obtained X-ray luminosity 
[26], within a broader beam that extends beyond the nucleus, is in every case already less than 10 -5 (in 
Eddington units). The four QRs that appear to be correlated with UHECRs (> 40 EeV) are identified 
here by their offset angles relative to the measured arrival direction of the primary particle initiating the 
associated air shower; this offset is never greater than the estimated experimental error in determining this 
arrival direction. For the worst case, NGC 4589, the associated Haverah Park UHECR direction measured 
has an angular uncertainty radius of 5.6 degrees. 

Along with the six published values cited in Table 4, we estimate the black hole mass (Mbh) f° r a U 
twelve QRs from the velocity dispersion of stars []27j within the bulge. The galactic bulge mass (Mqb) is 
derived from its stellar luminosity using the mass- luminosity formula (M/M©) / (L/10 10 L Q ) @. For 
the SO galaxy NGC 2300, we assume a typical SO bulge-to-disk ratio of my (bulge) — my (total) = 0.60 [29]. 
We note that the bulge mass loss rate in giant elliptical galaxies is ~ [Mqb / (W 12 M & )] M /year |K]] and 
assume that the rate of accretion onto the central black hole is an order of magnitude less, as typically found 
for the Bondi rate fr3l|| ; based on this, we estimate that the accretion rate is dM/dt ~ O.IM12 (M Q /year), 
where M12 = [(Mgb)/(10 12 M q )]. In Eddington units, this rate is then 

m = c 2 (dM/dt)/L Edd « 0.45(M 12 /M 8 ), (2) 

where M 8 ee [(M B h)/(1O 8 M )], and L Edd = 1.3 x 10 46 M 8 ergs/s is the Eddington luminosity limit. 



From equations 2 — 4 in the last paper of Ref. | |10[ |, we obtain that, for losses dominated by curvature 
radiation, the maximum proton energy expected via the dynamo action considered is 

£max = 77 (dM/dt) 1 / 8 (M 8 ) 1/4 EeV « 58 (M 12 ) 1/8 (M 8 ) 1/4 EeV. (3) 

We note that the values of E mSuX for all 12 galaxies in Table 4 lie above 40 EeV, the lower limit character- 
izing the AGASA sample considered. 

A lower limit to the radiation length (A m j n ) for proton energy loss associated with photo-pion 
production is estimated by considering the population of target photons within the source region 
[i?(source radius) > 2GM/c 2 ] at radio frequencies v > 360(7/10 11 )~ 1 GHz p0| , where 7 is the Lorentz 
factor given by 7 = (£' max )/(938MeV). The radio estimates required are extrapolated from data at lower 
frequencies fp2|. Those QRs with A m i n > R are expected to successfully accelerate protons up to ~ E max . 



Apart from their apparent correlation with UHECRs, there is no obvious systematic difference between 
the first four QR candidates listed at the top of Table 4 and the eight remaining ones at the bottom. 
However, it's important to note that we do not as yet know the spin states of the supermassive black holes 
associated with the nuclei of these galaxies; this key parameter might differ substantially among them. 
In each instance, the present state of spin depends on the specific prior history involved (e.g., accretion 
evolution, merger interactions, and earlier activity). Since the galactic nuclei considered here are X-ray 
dark, as are most (see first Ref. in |[l)), they are not viable candidates for black hole spin determination by 
means of the accretion disk iron K-line x-ray florescence that appears to be so promising a spectroscopic 
tool for active galactic nuclei, particularly Seyferts We emphasize that a black hole state of near 

maximal spin is a necessary condition for the realization of a compact black hole dynamo of the sort 
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Table 4: Physical attributes of quasar remnant candidates. Col. (1), catalog number. Col. (2), morpholog- 
ical type. Col. (3), distance, D, corresponding to the corrected redshift in Table 1 (for Hq = 75). Col. (4), 
galactic bulge mass, -Mgb, in units of 1O 12 M0. Col. (5), black hole mass, Mbh> in units of 1O 8 M0. Col. 
(6), accretion rate, m, in Eddington units. Col. (7), maximum proton energy, -E max - Col. (8), minimum 
radiation length, A m i n , relative to source size, R. Col. (9), offset angle between the putative QR and the 
most nearly aligned UHECR observed (> 40 EeV). 
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Note — Black hole masses in the table are normally calculated from stellar velocity dispersions. Those 
black hole masses taken from the literature are shown in parentheses, as are the quantities derived from 
them. A f symbol stands for masses taken from van der Marel 1999, whereas a X symbol, for masses taken 
from Kormendy & Gebhardt 2001 ||35|| . These literature masses are corrected for our assumed distances. 



envisaged for accelerating UHECRs [10|. In this sense, if confirmed, the correlation of UHECRs and QRs 



might well signal the introduction of a new means for identifying those nearby isolated non-active galactic 



nuclei that harbor highly spun- up black holes [34]. 



4 Concluding remarks 

As recently suggested, quasar remnants are plausible sites of black- hole based compact dynamos that could 
be capable of accelerating protons up to ultra-relativistic energies. We have found that nearby quasar 
remnant candidates present an above-random positional correlation with the sample of UHECRs. The 
correlation appears on closer angular scales than those expected when taking into account the deflection 
caused by typically assumed intergalactic or Galactic magnetic fields. Possible scenarios producing this 
effect were discussed; if real, the plausible fine structure of the Galactic field may ultimately provide the 
basis for the most natural explanation. In order to substantiate and further investigate the apparent 
correlation reported here between QR candidates and CR arrival directions, we need a large, reliable 
sample of the most energetic UHECRs, those of the very highest magnetic rigidity. It is hard to claim 
a definitive correlation with few objects and CRs constituting the samples. Future experiments, such as 
the Pierre Auger observatory [|36[ , EUSO [37] and the NASA space-borne OWL (Orbiting Wide-angle 



Light-collectors) mission []38|, should vastly increase the availability of UHECRs > 100 EeV. As we have 
already noted, however, CR sources in rich galactic clusters (those with pervasive micro-Gauss fields) are 
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not well suited for our present kind of correlative investigation, even though the nearest QR candidates of 
interest reside in such systems. We suggest that QR candidates located within the relatively nearby Virgo 
and Fornax clusters [ 10 1 might be best studied by means of the TeV curvature radiation expected from 



the putative compact black-hole dynamos associated with these objects (see second paper in [ 10 J ) . This 
work suggests that QRs should also be taken into account when analyzing coincidences in the forthcoming 
Auger Observatory. 
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